Fermionic stabilization and density-wave ground state of a polar condensate 
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We examine the stability of a trapped dipolar condensate mixed with a single-component fermion 
gas at T = 0. Whereas pure dipolar condensates with small s-wave interaction are unstable even 
for small dipole-dipole interaction strength, we find that the admixture of fermions can significantly 
stabilize them, depending on the strength of the boson-fermion interaction. Within the stable regime 
we find a region where a ground state is characterized by a density wave along the soft trap direction. 
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The recent demonstration of a condensate of chromium 
atoms [H opens up the study of quantum-degenerate 
gases that interact via long range, anisotropic magnetic 
dipole interactions. In a parallel development, it can be 
expected that quantum degenerate samples of heteronu- 
clear polar molecules will soon be available through the 
use of Feshbach resonances [l, Q , photoassociation [1, Q , 
or a combination of the two approaches, other possible 
routes including buffer-gas cooling collisional beam 
cooling Q, reactive scattering Q and stark deceleration 
alone or in combination. When in their vibrational 
ground state these molecules interact primarily via the 
electric dipole interaction, and may therefore also form 
dipole-dominated condensate. 

As a result of the dipole-dipole interaction, a number 
of novel phenomena have been predicted to occur in low- 
density quantum-degenerate atomic and molecular sys- 
tems, both in conventional traps and in optical lattices. 
The existence of a variety of quantum phases including a 
"supersolid" phase for dipolar bosons in optical lattices 
have been predicted in Rcf. [l^. A supersolid phase was 
also predicted for hardcore bosons in two-dimensional tri- 
angular lattices [ll|, , in Kagome lattices , and in 
an extended multi-band Bose- Hubbard Hamiltonian Il4i . 
although recent quantum Monte Carlo simulations [15[ 
did not find any supersolid phase for bosons in a Kagome 
lattice. The authors of Ref. studied Bose- Fermi mix- 
tures in two-dimensional square lattices, and found a 
bosonic supersolid transition induced by a modulation 
of the fermionic density resulting from a nesting effect. 
However, the situation is much different in the absence 
of a lattice structure: the only stable state of dipolar 
condensates in a pancake-shaped trap has a Gaussian- 
like density profile [13, [H, and states with a pe- 
riodic density modulation are always unstable [20j . al- 
though Ref. [2l| has shown that modulating periodically 
the strength of a laser-induced dipole-dipole interaction 
can result in stable density-modulated condensates. 

In this note we consider a mixture of dipolar bosons 
and non-interacting single-component fermions confined 
to a cigar-shaped trap. Such a mixture might be realized 
using isotopes of chromium, or during the production of 
heteronuclear molecules via either a Feshbach resonance 
or the photoassociation of two different fermionic species 



of atoms. Using linear response theory to determine the 
fermion-induced interaction between bosons we find that 
it significantly stabilizes the dipolar bosonic condensate. 
Wc then identify a region in parameter space where the 
stable ground state displays a density-modulated struc- 
ture along the long axis of the trap. 

We assume for simplicity that both the dipolar bosons 
and the fermions are trapped in the transverse direction 
by a tight harmonic potential of frequency lu± and in the 
longitudinal direction by a much softer harmonic poten- 
tial of frequency lOz ■ The dipoles are taken to be polarized 
by an electric or magnetic field in a direction y perpen- 
dicular to the long axis of the trap. The dipole-dipole 
interaction between two bosonic particles separated by 
a distance r is given by Vdd('') = ffdd (l — •iy^/r^) /r^, 
where gdd is the dipole-dipole interaction strength. In 
the mean-field approximation, the energy functional for 
the order parameter (p{r) of the dipolar condensate is 
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where Hq = -ti^V^ /{2m) + mu^ {x^ + + A^z^) /2 is 
the sum of the kinetic energy and the trapping poten- 
tial, m is the mass and the total number of bosonic 
particles, and A = ujz/u!_\_. the second term in E is the 
boson-boson s-wave interaction of strength g, and the 
third term describes the nonlocal dipole-dipole interac- 
tion between bosons. Note that that term also contains 
a short-range contact contribution [l^, [H, [2^ . The last 
term -Bind accounts for the fermion-induced interaction 
14nd(fe) between bosons. It is given by 
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Vind{k)n{k)n{~k)d^k, 
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where n{k) is the momentum-space bosonic density. 

In deriving Eq. ([2]) we have assumed a contact boson- 
fermion interaction of strength gbf, with a correspond- 
ing interaction energy of the form (/bf J nf{k)n(~k)d^k, 
nf{k) being the fermion density. The linear response of 
the fermions to a bosonic density fluctuation n{k) can 
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be expressed as nf{k) = Vind{k)n{k), and Eq. ([2]) is ob- 
tained by substituting that expression back into to the 
boson-fermion interaction energy. The explicit form of 
the induced potential is V^nd(fc) = Q hfXfjk), where Xf 
is the density response function [25| related to the dy- 
namical structure factor S{k,uj), the probability of ex- 
citing particle-hole pairs with momentum k out of the 
Fermi sea, by x/(fc) = -2 duj'[S{k,uj')/uj']. For a 
non-interacting single-component Fermi system we have 
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Here p = \p\, kf is the Fermi momentum, and the excita- 
tion energy is ujp^ = pk cos6'/m/ + fc^/ {2m f), 9 being the 
relative angle between p and fc, and to/ the mass of the 
fermions. This expression assumes that the fermions are 
locally free, so that there is a Fermi sphere in momentum 
space. 

Using the form of S{k,ijj) from Ref. [1^ we find 
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where ly is the three-dimensional density of states at 
the Fermi surface. Here |..| is used to denote abso- 
lute value. Substituting Eq. ([4]) into the expression 
Vindik) = ghfXfi^) carrying out a series expansion 
to lowest order in the non-local terms gives 
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(5) 



The negative contribution to this potential describes an 
effective attractive interaction between bosons in the long 
wavelength limit, while its positive part corresponds to a 
nonlocal repulsive interaction. 

With the form ^ of the bosonic energy functional 
and the potential Vind at hand, we now proceed to de- 
termine the stability of the dipolar condensate, using a 
variational wave function in the parameter space of the 
fermion-induced interaction and dipolar strength. For 
this purpose, the condensate order parameter (f>{r) is as- 
sumed to factorize as (/)(r) = </>_L(a;, y)(f>\\ (z) with the nor- 
malization J (fr\(j){r)\'^ = 1. The transverse wave func- 
tion is (j)±{x,y) = exp[— (x^ -I- y^) / {2d'^)]/ VtuP, d being 
a variational parameter, while the longitudinal compo- 
nent is (jiwiz) = exp[— z^/(2d^)]/v7r'72rf^, dz being like- 
wise a variational parameter. 

Substituting this Gaussian ansatz and its Fourier 
transform into Eqs. H]) and the energy of the con- 
densate becomes 
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FIG. 1: Stability diagram of a dipolar condensate in ((/3d, 5ind) 
space for Qs = 0.2 and A = 0.05. In region I the stable ground 
state has a Gaussian-like density distribution and in region II 
it is characterized by a density modulation, i.e., adw > 0. 
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where 77 = d/d^, kf ~ kfd±, d± = y^h/muj± is the 
transverse oscillator length, and 
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The first term in Eq. ^ denotes the short range part 
of Vddi and the second term results from the difference 
between the boson-boson s-wave interaction and the at- 
tractive part of the interaction in Vind{k)- To calculate 
the fermion-induced interaction energy Ei^d we have as- 
sumed that 2k f > I. In Eq. dSl) we have introduced 
the effective three-dimensional dipole-dipole interaction 
ff3d = N gddfn/ {\/2nh? d±) , the effective contact potential 
gs = l/3+(g— (7^fi^)/(27rgdd), and the induced interaction 
5ind = glfiy/{l2TTgddkj)- 

The energy functional ^ is minimized with respect 
to d/d± and j] — d/ dz for various combinations of the 
parameters g^d, 5ind- The results of these simulations are 
summarized in Fig. [TJ All numerical results presented in 
the following arc for A = 0.05, gs = 0.2, and kf = 2 
(a value typical of the situation for a fermionic gas of 
10^ atoms in a a transverse trap length of d± ^ l/im 
and A ~ 0.05). We found that the condensate typically 
collapses for 77 > 1 or d ^ 0, while the stable regime 
is characterized by 77 < 0.1 and d ^ d±. Region I is 
characterized by stable ground states with a Gaussian- 
like density profile, while in Region II the ground state 
exhibits density modulations that are further discussed 
later on. 

In the absence of fermion-induced effective boson- 
boson interaction g-^d = 0, the dipolar condensate is 
found to be unstable to collapse even for very small dipo- 
lar strength g^d- For ^ind 7^ 0, in contrast, the conden- 
sate becomes stable for a range of values of gsd that in- 
creases with increasing gind- This can be understood from 
the form of the non-local part of the induced interac- 
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FIG. 2: Effective one-dimensional potential, Eq. (jS)), for sev- 
eral values of the fermion-induced interaction (?ind- Here we 
assumed that d ~ dx- The zero potential level is shown as a 
horizontal dashed line. 



FIG. 3: Qualitative diagram of the Bogoliubov spectrum for 
the effective longitudinal potential from Eq. (jSj for several 
values of gsd and an arbitrary finite value of gind. 



trum 



tion, Eq. ([5]). Since the induced interaction Vi^d is repul- 
sive for high momenta, see Fig. 3, the non-local repulsive 
part balances the attractive non-local part of the dipole- 
dipole interaction with increasing ^ind, which results in 
increased stability of the condensate. We confirmed that 
the qualitative features of the stability diagram remain 
unchanged for 0.05 ^ 5s 0.3. The system is also found 
to become stable for higher dipolar strengths gsd with 
increasing g^. 

We now restrict our considerations to the case of a 
cylindrical trap, and discuss the appearance of a density- 
modulated ground state in region II of the stability dia- 
gram. We proceed by first demonstrating the existence 
of a roton minimum in the Bogoliubov excitation spec- 
trum of a homogeneous condensate along the trap axis, 
and show that by increasing the effective dipolar strength 
(73d this minimum can touch the zero-energy axis. The 
roton-like spectrum follows from the attractive nature of 
dipolar interaction for high momenta [H, [131 and 

touching the zero energy indicates an instability toward 
a nonuniform ground state [28j . 

Integrating Eqs. ^ and ^ in the transverse direction 
yields the effective longitudinal potential 



Vcs{kz) = 9s- 




9ind^kl (8) 



where kz = fczd is a scaled wave vector, kz being the lon- 
gitudinal momentum, and Ei{x) = dx{e~^ / x) is the 
exponential integral. In obtaining Eq. ([8|) we have as- 
sumed that kz < kf which justifies keeping only the con- 
tribution to the induced potential ([5]) with k < 2k f. Fig- 
ure [2] shows that as gind is increased the one-dimensional 
potential becomes less attractive for higher momenta, 
and also that the non-local part of Vcf[{kz) becomes re- 
pulsive again for sufficiently strong g-md- 

Assuming a homogeneous condensate of length L along 
the long axis of the trap then yields the Bogoliubov spec- 
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with the appearance of a roton minimum for large enough 
effective dipolar strength, see Fig. [3) By increasing g^d 
further, the roton minima can touch the zero-energy axis. 
This suggests that at this point the structure of the con- 
densate ground state can undergo a transition, with the 
appearance of a density modulation along the long axis 
of the trap. 

To verify this possibility we introduce a new longitu- 
dinal variational wave function that includes a density- 
wave component of amplitude Odw 



,dw(^) = 0||(z) 



a + fldw cos 




with the normalization constraint 



(9) 



J2 = 1. The 



excess energy of this density-modulated state, as com- 
pared to a broad Gaussian order parameter, can be ex- 
pressed as 



-- 2mtP 
o r 2 



(10) 



'^ + ^93d^\ 2aVeff(fco) + ^14ff(2fco) 



Limiting ourselves to the stable region of Fig. 1, we 
first choose the values of d and dz that minimize the en- 
ergy ([HI), and then minimize e(fco, ftdw) with respect to the 
variational parameters fcg and Odw Density- wave ground 
state is characterized by min {e} < and Odw 0. 
The region of density-wave ground state determined in 
this fashion corresponds to region II of Fig. [1] For 
.9ind > 0.023, a transition from a Gaussian-like state to 
a density-wave state is observed when the dipole-dipole 
interaction g^d is increased. This transition is smooth 
in the sense that Odw increases gradually from zero as 
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FIG. 4: Longitudinal density profile of the condensate as a 
function of increasing dipolar strength g^^ for gind = 0.058 
inside region II of Fig. 1. The longitudinal displacement z is 
scaled to the transverse oscillator length d±. (a) g^d = 7.9, 
where the system just enters region II, with excess energy 
e = -1.0 x'lO"^ (b) p3d = 12, e = -3.5 x 10"^ and (c) 
g3d = 20, e = -2.0 X 10"^ 



we increase the dipolar strength g^^ beyond the critical 
value. Figure m illustrates the gradual emergence of den- 
sity waves for a fixed induced strength gind = 0.058. We 
find that by gradually increasing (73d to a modest value 
the density profile acquires periodic zeros along the lon- 
gitudinal direction. 

One possible candidate to observe the predicted den- 
sity wave region is a mixture of bosonic^^Cr and fermionic 
^^Cr atoms. Bosonic Chromium has a magnetic dipole 
moment of and a s-wave scattering length ^ IOOoq 
[29I . [30| . The trap considered in this paper has a trans- 
verse frequency u}± ~ 200Hz and aspect ratio A = .05. 
For this geometry the density wave state in region II 
of Fig. [1] can be reached for a boson-fcrmion scattering 



length of 500ao. for bosonic and fermionic atoms num- 
bers of 10"' and 10'^, respectively. 

In summary, we have analyzed the stability of ultra- 
cold dipolar bosons mixed with non-interacting fermions. 
A central result of our analysis is that the fermions help 
stabilize the dipolar condensate, a consequence of the 
fact that the non-local induced interaction is repulsive 
in the limit of moderate wavelengths. In the stable re- 
gion we found a transition in the shape of the condensate 
ground-state from a Gaussian-like profile to a modulated 
density profile along the axis of the trap. Our analy- 
sis is mean-field, implying that the system is assumed 
to be condensed with phonon-like low energy excitations 
along the longitudinal direction. The density modula- 
tion emerges as a result of the additional breaking of 
the translational symmetry of the condensate and can be 
called a superfluid density wave. 

Future work will discuss the effect of the strength of 
the contact interaction and of the Fermi momentum on 
the stability of the condensate and the density wave, with 
possible extensions to pancake geometries and to rotat- 
ing systems. A number of exotic states are likely to be 
found in that regime. An extension to finite temperatures 
will allow us to investigate possible classical phase tran- 
sitions between superfluid and various kinds of density- 
wave states [Slj. 

This work is supported in part by the US Office of 
Naval Research, by the National Science Foundation, and 
by the US Army Research Office. 
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